Antiferromagnetic Ising model in scale-free networks 
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The antiferromagnetic Ising model in uncorrelated scale-free networks has been studied by means 
of Monte Carlo simulations. These networks are characterized by a connectivity (or degree) distribu- 
tion P{k) ~ The disorder present in these complex networks frustrates the antiferromagnetic 
spin ordering, giving rise to a spin-glass (SG) phase at low temperature. The paramagnetic-SG 
transition temperature Tc has been studied as a function of the parameter 7 and the minimum 
degree present in the networks. Tc is found to increase when the exponent 7 is reduced, in line with 
a larger frustration caused by the presence of nodes with higher degree. 
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I. INTRODUCTION 

Several types of natural and artificial systems have a 
network structure, where nodes represent typical system 
units and edges play the role of interactions between con- 
nected pairs of units. This kind of description of complex 
systems as networks or graphs has attracted much inter- 
est in recent years. Thus, complex networks have been 
used to model various types of real-life systems (biologi- 
cal, social, economic, technological), and to study several 
processes taking place on them [H, S Hi ■ Some mod- 
els of networks have been designed to explain empirical 
data in various fields. This is the case of the so-called 
small- world and scale- free networks which incor- 
porate different aspects of real systems. 

In scale-free (SF) networks the degree distribution 
P(fc), where k is the number of links connected to a node, 
has a power-law decay P{k) ^ k'~^ [1, [§]. This kind 
of networks have been found in the internet llOll , in the 
world-wide we b for protein interactions and in 
social systems [Ij]- The origin of such power-law degree 
distributions was addressed by Barabasi and Albert [tI, 
who found that two ingredients are sufficient to explain 
the scale-free character of networks: growth and prefer- 
ential attachment. They concluded that the combination 
of these criteria yields non-equilibrium SF networks with 
an exponent 7 = 3. One can also consider equilibrium 
SF networks, defined as statistical ensembles of random 
networks with a given degree distribution P{k) [1, [13, 
for which it is possible to analyze various properties as a 
function of the exponent 7. 

Cooperative phenomena in complex networks display 
unusual characteristics due to their peculiar topology 
[H [H, [13, [ll E, H^l- In particular, the ferromagnetic 
(FM) Ising model has been thoroughly studied in scale- 
free networks [ll], |23, [H, [IJl, and its critical behavior 
was found to be dependent on the exponent 7. For finite 
(/c^), a ferromagnetic to paramagnetic transition occurs 
at a finite temperature T^. However, when (fc^) diverges 
(as happens for 7 < 3), the spin system remains in an 
ordered FM phase at any temperature, so that no phase 
transition appears in the thermodynamic limit. 



Here we study the antiferromagnetic (AFM) Ising 
model in equilibrium (uncorrelated) scale-free networks 
with several values of the exponent 7. This model con- 
tains the two basic ingredients necessary to produce 
a spin-glass (SG) phase at low temperature: disorder 
and frustration. In some spin-glass models, such as the 
Sherrington-Kirkpatrick model, all spins are assumed to 
be mutually connected [2^, [2y|, whereas in others ran- 
dom graphs with finite (low) connectivity are considered 
[13, [3, IM l30|. For the AFM Ising model on scale- 
free networks, we expect to find features intermediate 
between these two cases. 

Spin glasses on complex networks have been studied in 
recent years by usirig several techniques, such as trans- 
fer matrix analysis [3l| , replica symmetry breaking [3^ , 
defect- wall calculations [33| , and an effective field theory 
[33 . In this paper, we employ Monte Carlo (MC) sim- 
ulations to study the paramagnetic to spin-glass phase 
transition appearing in scale-free networks. In this line, 
MC simulations have been carried out earlier to analyze 
spin-glass phases appearing for the AFM Ising model in 
Barabasi- Albert scale-free networks [35|, as well as in 
small- world networks [3^ . 

The paper is organized as follows. In Sec. II we de- 
scribe the networks and the computational method used 
here. In Sec. Ill we present results for the heat capacity, 
energy, and spin correlation, as derived from MC simu- 
lations. In Sec. IV we characterize the spin-glass phase 
through the overlap parameter and transition tempera- 
ture. The paper closes with the conclusions in Sec. V. 



II. MODEL AND METHOD 

We consider SF networks with degree distribution 
P{k) k~'^ . They are characterized, apart from the ex- 
ponent 7 and the system size N , by the minimum degree 
/co- We assume that P(k) = for fc < fco. Our networks 
are uncorrelated, in the sense that degrees of nearest 
neighbors are statistically independent. This means that 
the distribution P(fc, k') of degrees of nearest-neighbor 
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nodes fulfills the relation [J 

kk' 



P{k,k') 



{ky 



p{k)p{k') 



(1) 



Alternatively, one can use a correlation coefRcient r de- 
fined as 



{kk') - {k){k') 



(2) 



where the averages are taken over all links and of, is the 
variance of the degree distribution. This coefficient r is 
zero for uncorrelated networks. 

For the numerical simulations we have generated net- 
works with several values of 7, /cqj and N . To generate 
a network, we first define the number of nodes with 
degree fc, according to the distribution P{k), which can 
be conveniently done by using the so-called transforma- 
tion method [33] ■ Then, we ascribe a degree to each node 
according to the set {Nk}, and finally connect at random 
ends of links (giving a total of i = kNk/2 connec- 
tions), with the conditions: (i) no two nodes can have 
more than one bond connecting them, and (ii) no node 
can be connected by a link to itself. We have checked that 
networks generated in this way are uncorrelated, i.e. they 
fulfill Eq. ([1]), and r — 0. The networks considered here 
contain a single component, i.e. any node in a network 
can be reached from any other node by traveling through 
a finite number of links. 

Given a network with a particular set of links, we con- 
sider an Ising model with the Hamiltonian: 



H — ^ JijSiSj , 



(3) 



where Si — ±1 {i = l,...,iV), and the coupling matrix 
Jij is given by 



J,. 



_ / J(> 0), if I and j are connected, 
0, otherwise. 



(4) 



This means that each edge in the network represents an 
AFM interaction between spins on the two linked nodes. 
Note that, contrary to the usually studied models for 
spin glasses in which both FM and AFM couplings are 
present, in our model all couplings are antiferromagnetic 
(similarly to Refs. [Hill). 

For a given network, we carried out Monte Carlo sim- 
ulations at several temperatures, sampling the spin con- 
figuration space by the Metropolis update algorithm [39*1 , 
and using a simulated annealing procedure. Several vari- 
ables characterizing the spin system have been calculated 
and averaged for different values of the considered pa- 
rameters. For each set of parameters (7, fco, N), 1000 
network realizations were generated, and the largest net- 
works included 16000 nodes. In the sequel, we will use 
the notation (...) to indicate a thermal average for a net- 
work, and [...] for an average over networks with a given 
parameter set. 




FIG. 1: Heat capacity per node, c„, vs temperature for scale- 
free networks with A'^ — 8000 nodes and minimum degree 
fco = 5. The plotted curves correspond to different values 
of the exponent 7: 10 (solid line), 5 (dashed line), 4 (dotted 
line) , and 3 (dashed-dotted line) . 



III. THERMODYNAMIC VARIABLES 

We present first results for the heat capacity per site, 
Cv, as a function of temperature for several values of the 
exponent 7. c„ has been derived from the energy fluctua- 
tions AE at a given temperature, by using the expression 



(5) 



where (AE)'^ = {E"^) - (E)'^. We have checked that the 
results obtained in this way coincide within statistical 
noise with those derived by calculating the heat capacity 
as the energy derivative [d{E) / dT] / N . Note that we take 
the Boltzmann constant fc^ = 1. 

The temperature dependence of Ci, is displayed in Fig.[T] 
for scale-free networks with various values of 7 between 3 
and 10. The data shown correspond to networks includ- 
ing 8000 nodes. For increasing 7, we observe that the 
maximum of Cy shifts to lower T, and the peak becomes 
narrower. This narrowing is in line with that observed 
earlier for the heat capacity in the AFM Ising model on 
small- world networks, when the disorder is reduced [3^ . 
In our case of scale-free networks, larger values of the 
exponent 7 correspond to networks with a higher homo- 
geneity in the node connectivities (less dispersion in the 
degree distribution), causing a narrower peak in the heat 
capacity, as shown in Fig. 1. The peak shift to lower tem- 
peratures suggests a transition from a paramagnetic to a 
SG phase at a temperature that decreases as the expo- 
nent 7 rises. For increasing 7, one reduces the presence 
of nodes with a large degree, which in turn reduces the 
degree of frustration in the spin distribution (see below). 
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FIG. 2: Minimum energy per link derived from our simula- 
tions for the AFM Ising model on scale- free networks, plotted 
as a function of the inverse exponent 7~^. Symbols represent 
energy values obtained in the large-network limit N ^ oo 
for several values of the exponent 7, and minimum degree: 
ko = 10 (triangles); fco = 5 (circles); fco = 3 (squares). Sym- 
bols at — correspond to regular random networks with 
a constant degree. Lines are guides to the eye. 



AFM ordering on the considered random networks 
with power-law distribution of degrees is frustrated by 
the disorder in the hnk configuration, and in particular 
by the presence of loops with an odd number of nodes. 
The degree of frustration can be quantified by looking 
at the low-temperature energy of the system, which will 
be higher for larger frustration. Given the parameters 7 
and ko defining the scale-free networks, we obtain a value 
for the minimum energy by extrapolating to infinite size 
the minimum energy reached in our simulations of finite- 
N networks. This extrapolation has been performed by 
assuming a dependence of the energy on network size of 
the form: 



iV2/3 ' 



(6) 



where em{N) is the energy per link for size N, Cm is its 
limit for N 00, and A is a fit parameter. This kind of 
dependence of the low-temperature energy in spin-glass 
systems was proposed by Boettcher [s^, and has been 
found to be followed by the results of our calculations 
using the simulated-annealing method. We have checked 
that our method to obtain a minimum energy for the 
AFM Ising model on complex networks gives similar re- 
sults to those found by using extremal optimization [soj ■ 
In particular, for spin glasses on random graphs with a 
Poisson distribution of connectivities, we found results 
very close to those obtained by using this technique [s^ . 
In any case, the energy em found here for each parameter 
set {ko, 7) will be an upper limit for the lowest energy of 
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FIG. 3: Spin correlation function ^ vs distance for scale-free 
networks with fco = 3 and several values of the exponent 7, at 
temperature T — J. Symbols correspond to different values 
of 7: 10 (squares), 5 (circles), 4 (triangles), and 3 (diamonds). 
These results were derived from simulations for SF networks 
including A*" = 8000 nodes. 



the system. 

In Fig. [2] we show results for the minimum energy per 
link Cm found for three values of fco and several values of 
the exponent 7. For a given fco one observes an increase 
the exponent 7 is reduced. This indicates that 
the presence of nodes with large degree (hubs), which 
is favored for small 7, plays in our context the role of 
increasing the frustration in the spin arrangement. For 
a given 7, one observes also in Fig. [2] an increase in the 
energy Cm for rising fco. This shows that an increase in 
the minimum connectivity (or in the average degree (fc)), 
causes also a higher frustration in the AFM ordering. 

A quantification of the short-range order present in 
the spin system on scale-free networks can be obtained 
by calculating the spin correlation 



e(r) = [{S,S,)r 



(7) 



where the subscript r indicates that the average is taken 
for the ensemble of pairs of sites at distance r. Note 
that the dimcnsionless distance r refers to the minimum 
number of links between two nodes, also called in the lit- 
erature chemical or topological distance. The correlation 
^(r) is shown in Fig.[3]for several values of the exponent 7 
at a temperature T — J. This temperature is below the 
critical temperature Tc of the paramagnetic-SG transi- 
tion for all values of 7 (see below). As expected, ^(r) de- 
creases faster for smaller 7, due to the presence of nodes 
with large degree, and consequently a larger frustration 
of the AFM ordering, as discussed above in connection 
with the minimum energy per link. 

To obtain more direct insight into the reduction of 
^(r) with the distance, we display in Fig. [5] \£,{r)\ on a 
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FIG. 4: Absolute value of the spin correlation function ^ 
vs distance in a semilogarithmic plot. Shown is ln|^(r)| for 
scale-free networks with minimum degree fco = 3 and several 
values of the exponent 7 at T = J. Symbols correspond to 
different values of 7. From top to bottom: 7 = 00 (regular 
networks, open squares), 10 (circles), 6 (triangles), 5 (open 
diamonds), 4 (filled diamonds), and 3 (filled squares). Error 
bars are on the order of the symbol size. These results were 
derived from Monte Carlo simulations for networks including 
iV = 8000 nodes. 



semilogarithmic plot for various 7 values. In general, 
after a short transient for small r, one obtains an ex- 
ponential decrease in the spin correlation with distance, 
as |^(r)| ~ e"'"', a being a parameter that depends on 
temperature as well as on the parameters defining the 
networks (fco and 7). For the results shown in Fig. |4l we 
find a parameter a that decreases from 1.01(3) to 0.386(5) 
when increasing 7 from 3 to 10, and reaches the limit a 
— 0.321(4) for regular random networks with constant 
degree k ^ ko = 3. We note that here the limit 7 ^ cxd 
correspond to networks (called regular [Soj), in which all 
nodes have the same degree fcg, and consequently do not 
include any hub with high degree. 
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FIG. 5: Distribution of the overlap parameter q for scale- 
free networks with minimum degree fco = 5, as derived from 
MC simulations, (a) Networks with exponent 7 = 3 at three 
temperatures: T/J — 4, 3, and 2; (b) networks with three 
different exponents 7 at temperature T = 2J. 




IV. SPIN-GLASS BEHAVIOR 

A. Overlap parameter 

In the study of spin glasses, it is usual to consider two 
copies of the same network, with a given realization of the 
disorder. Then, one considers a spin system on each net- 
work, both with different initial values of the spins, and 
follows their evolution with different random numbers for 
generating the spin flips [U, \^ . A particularly relevant 
parameter is the overlap q between the two copies, de- 



fined as 

i 

where the superscripts (1) and (2) indicate the copies. 
This parameter q is defined in the interval [—1,1], and the 
extreme values 1 and -1 correspond to pairs of networks 
with the same spin configuration (apart from a trivial 
overall flip in the -1 case). 

We have calculated the overlap parameter q for scale- 
free networks with various exponents 7, and derived the 
probability distribution Q{q) from Monte Carlo Simula- 
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tions. Results for Q{q) are presented in Fig. [5l In par- 
ticular, in Fig. 5(a) we display the distribution of the 
overlap parameter for networks with 7 = 3 at several 
temperatures. At high temperatures (T ^ J), the distri- 
bution Q{q) shows a single peak centered at q — 0, which 
is characteristic of a paramagnetic state. This peak has, 
however, a finite width, which results to be a finite-size ef- 
fect. It should collapse to a Dirac delta function at g = 
in the limit N 00. When the temperature is lowered, 
Q{q) broadens around q = 0, due to the appearance of 
an increasing number of frustrated links. At still lower 
temperatures, frustration is more apparent, and "freez- 
ing" of the spins causes the appearance of two peaks in 
Q{q), symmetric respect to q = 0, and characteristic of 
spin-glasses [ssl, fS, l43i| . Such a distribution Q{q) is 
associated to the break of ergodicity occurring in the spin 
system at low temperatures. 

In Fig. 5(b) we show the distribution Q{q) for three 
values of the exponent 7 at a fixed temperature T = 2 J. 
The effect of decreasing 7 for a given T is similar to that 
shown in Fig. 5(a) for lowering the temperature for a 
given 7, in the sense that in both cases one passes from 
a high-temperature paramagnetic phase to a spin-glass 
with broken ergodicity. From the results displayed in 
Fig. 5(b), along with those presented in Sect. Ill (spe- 
cially Fig. 1 for the heat capacity and Fig. 2 for the mini- 
mum energy per link) , we expect that the freezing of the 
spins in the SG phase occurs at lower T for larger 7. In 
other words, one expects that the transition temperature 
from paramagnetic to SG will decrease for rising 7. 



B. Transition temperature 

The overlap parameter q can be further employed to 
obtain precise values of the paramagnetic-SG transition 
temperature. To this end, one can use the fourth-order 
Binder cumulant [39, 531 



5iv(T) 



3- 



N 



(9) 



N 



which is restricted to the interval [0,1]. On one side, this 
parameter vanishes at high temperatures, as expected 
for a Gaussian distribution Q{q) in a paramagnetic state. 
On the other side, one has = 1 whenever the dis- 
tribution Q{q) vanishes everywhere except for \q\ = 1, 
which corresponds to the case of a single ground state, 
and could be reached at low temperatures. This case 
((/AT = 1) is clearly not expected here due to the onset of 
frustration, which gives rise to the spin-glass phase. 

In general gjv increases as temperature is lowered, and 
the transition temperature Tc can be obtained from the 
single crossing point for different network sizes N j36l.l39t. 
By employing this method, we have calculated Tc for 
scale-free networks with several values of the parameters 
ko and 7, and the results obtained are displayed in Fig.[6l 
In this figure, we have plotted Tc as a function of 
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FIG. 6: Transition temperature Tc from paramagnetic to 
SG phases, as a function of the inverse exponent 7~^, for 
three values of the minimum degree feg. Error bars, when not 
shown, are on the order of the symbol size. Dashed lines are 
guides to the eye. 



and the limit 7 ^ cxd corresponds to regular random 
networks with k = ko. This limit gives a reasonable ex- 
trapolation of the Tc values obtained for 7 values up to 
7 = 10. For a given minimum degree ko, the transition 
temperature increases as 7 is lowered. Also, for a given 
value of 7, Tc grows for increasing feg, similarly to the 
case of the minimum energy e™ displayed in Fig. O We 
note that the transition temperature Tc derived from the 
Binder cumulant is in the order of the maximum (nega- 
tive) temperature derivative of the heat capacity Cu for 
finite networks, as can be seen by comparing results for 
/co = 5 in Figs. 1 and 6. 

The average degree (k) of scale-free networks can 
be estimated rather accurately by replacing the sum 
kP{k) by an integral. Thus, one finds for networks 
with P{k) ~ k and 7 > 2, that the average degree 
scales as 



{k) 



7 ■ 



1 . 



7 ■ 



(10) 



From this expression, it is clear that d{k)/dko > and 
d{k)/dj < 0. Since rising ko or lowering 7 cause an 
increase in the transition temperature, we observe that 
in fact a rise in (fc) is accompanied by a rise in Te- 
la this context, it is interesting to compare our re- 
sults for Tc in scale-free networks with those found for 
the paramagnetic-SG transition in other complex net- 
works. In Ref. i36| it was studied the AFM Ising model 
on small-world networks generated by rewiring links in 
a regular lattice @. It was found that the transition 
temperature decreases as the disorder (number of ran- 
dom connections) increases, in networks where the aver- 
age degree was kept constant ((fc) =4). In the limit of 
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large disorder, those networks approach random networks 
with a Poisson distribution of degrees, and the transition 
temperature for the AFM Ising model was found to be 
Tc = 1.71J. Going back to our random networks with 
power-law degree distribution, we have an average degree 
(fc) « 4 for fco = 3 and 7 = 5 [see Eq. ^j]. For these 
networks, we find Tc = 1.65 J, as shown in Fig. 6, a value 
close to that obtained for small-world-type networks in 
the large-disorder limit and (fc) = 4. 

All this could suggest that the transition temperature 
Tc is proportional to the average degree (fc) , irrespective 
of the details of the networks under consideration. This 
is, however, not the case, as can be inferred directly from 
the results displayed in Fig. 6. Looking for example at 
the data for fco = 3, we observe that going from 7 = 3 
to the limit 7 ^ cx), (fc) decreases by a factor of 2 (from 
6 to 3). On the other side, Tc is reduced by a factor 
K. 4. Since lowering 7 increases the inhomogeneity in the 
degree distribution, favoring the presence of nodes with 
degree much larger than (fc), we find that this inhomo- 
geneity helps to rise the transition temperature. 

Bartolozzi et al. [s^ studied the AFM Ising model 
on Barabasi-Albert scale-free networks, and found a 
paramagnetic-SG transition temperature Tc = 4.0(1)J. 
These nonequilibrium networks are characterized by an 
exponent 7 = 3, and those authors used the particular 
value of the minimum degree fco = 5. For these param- 
eters, we find for equilibrium networks a transition tem- 
perature Tc = 4. 8(2) J, a value somewhat higher than 
that obtained in Ref. ISSl 

We note that the error bar in the transition tempera- 
ture grows when 7 is reduced. In fact, the actual value 
of the cumulant ^at at the crossing point for different N 
values (network sizes) decreases as 7 is lowered, and is 
near zero for 7 3. This coincides with results shown 
for this cumulant in Ref. |3^ for Barabasi-Albert SF net- 
works, where the value oi gn at the crossing point was 
less than 0.01. This means that the signal-to- noise ratio 
in gjv becomes poor, and one has an increasing uncer- 
tainty in Tc- For networks with 7 < 3, we could not find 
a single crossing point for the cumulant corresponding 
to different network sizes, and a transition temperature 



cannot be given. We note that these 7 values correspond 
to SF networks with diverging (fc^). In this respect, it is 
known that the ferromagnetic Ising model in such net- 
works does not show a phase transition, and remains in 
an ordered FM phase at any temperature in the thermo- 
dynamic limit [U, [22, |2^ [231 . Something similar could 
happen for the AFM Ising model in these networks. This 
point remains as a challenge for future research. 



V. CONCLUSIONS 

The AFM Ising model in random networks with a 
power-law distribution of degrees gives rise to a spin- 
glass phase at low temperature. This is a consequence of 
the combination of disorder in the networks and frustra- 
tion caused by the presence of loops with odd number of 
links. The overlap parameter q gives us evidence of this 
frustration at low temperatures. 

The transition temperature Tc from the high- 
temperature paramagnetic phase to the spin-glass has 
been studied as a function of the minimum degree fco and 
the exponent 7 in the degree distribution. Tc is found to 
rise for increasing fco and for lowering 7. 

For a given fco , both the transition temperature and the 
minimum energy per link found from our simulations 
increase as the exponent 7 is lowered. This indicates that 
the degree of frustration in the spin configurations rises 
with the presence of nodes with large degree (hubs). The 
same conclusion can be reached by analyzing the spin 
correlation as a function of distance, which decays faster 
for smaller values of 7. 
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